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Phase diagram of half-doped perovskite manganites is studied within the extended double- 
exchange model. To demonstrate the role of orbital degrees of freedom both one- and two-orbital 
models are examined. A rich phase diagram is obtained in the mean-filed theory at zero temperature 
as a function of J (antiferromagnetic (AFM) superexchange interaction) and V (intersite Coulomb 
repulsion). For the one-orbital model a charge-ordered (CO) state appears at any value of V with 
different types of magnetic order which changes with increasing J from ferromagnetic (F) to AFM 
ones of the types A, C and G . The orbital degeneracy results in appearance of a new CE-type spin 
order that is favorable due to opening of the "dimerization" gap at the Fermi surface. In addition, 
the CO state appears only for V > Vc for F and CE states while C-type AFM state disappears and 
A-type AFM state is observed only at small values of 1^ as a charge disordered one. The relevance 
of our results to the experimental data are discussed. 

PACS numbers: 75.25.4-z, 75.30.Et, 75.30.Vn, 71.45.Lr 



I. INTRODUCTION 

Since early fiftiesEl the physics of manganites challenge 
our current understanding of transition-metal oxides, and 
define both theoretical and experimental research prob- 
lem that involves charge, spin, lattice and orbital degrees 
of freedom. Recently in a modern systematic experimen- 
tal studies a very rich phase diagram (see, for example, 
Ref. 0) depending on the doping concentration, tempera- 
ture and pressure was obtained in the doped manganese 
oxides with perovskite structure Ri_^Ba;Mn03 (where R 
is trivalent rare-earth and B is divalent alkaline ion, re- 
spectively). At different doping concentration a full vari- 
ety of magnetically ordered states such as antiferromag- 
netic (AFM) insulator, ferromagnetic (FM) metal and 
charge ordered (CO) insulator were observed. Many ef- 
forts have been made by theoreticians to understand it 
based on various models and|-approaches. Historically, 
double exchange (DE) modelO was the basic one. In 
this model t23-electrons are localized, whereas the Cg- 
electrons are mobile and use O p-orbitals as a bridge 
between Mn ions. The hopping of itinerant electrons to- 
gether with a very strong on-site Hund's coupling drives 
core spins to align parallel. Qualitatively DE model gave 
appropriate interpretation of the phase diagram at the 
doping range 0.2 < x < 0.5 where FM metallic behavior 
was observed. 

One of the subtle aspects of the perovskite manganites 
is the charge ordered states observed in almost all such 
compounds at half-doping. qt3 A direct evidence of the 
CO state in half-doped manganites has beeruprovided by 
the electron diffraction for Lao.5Cao.5MnO3.la Similar olt 
servations have also been reported for Pro.5Sro.5Mn03El 
and in Ndo.5Sro.5MnO3.Ll CO state is characterized by 
an alternating Mn and Mn ions arrangement in 
X ~ y plane with the charge stacking in z-direction. 
In CO state these systems show an insulating behav- 



ior with a very peculiar form of AFM spin order- 
ing. The observed magnetic structure is a CE-type and 
consists of quasi one-dimensional ferromagnetic zig-zag 
chains coupled antiferromagnetically. In addition, these 
systems show d^^^ -^^ / d3y2 _j.2 orbital ordering. An- 
other noteworthy observations were done by studying 
Pri_a;(Cai_j,Srj,)a;Mn03 crystals with controlled one- 
electron bandwidth. As already mentioned above at 
half-doping Pro.5Cao.5Mn03 has a CO CE-type insulat- 
ing state. However, by substitution Ca with Sr lead- 
ing to the increase of the carrier bandwidth, one in- 
duces the collapse of the CO insulating state, and the 
A-type metallic state witi d^2_y2 orbital ordering is re- 
alized in Pro.5Sro.5MnO3.El The coexistence of the A-type 
spin ordered and CE-type spin/charge ordeted states has 
been detected in the bilayer_LaSr2Mn207cl and three- 
dimensional Ndo.5Sro.5MnO3.El These results indicate the 
competition between the metallic A-type dx2-y2 orbital 
ordering and the insulating CE-type d^x2_^2 Id^y2_^2 or- 
bital ordering at half-doping and demonstrate the impor- 
tance of the magnetic, charge and orbital order coupling 
in these compounds. 

In recent publications it was shown that the DE 
anisotropy resulted from the orbital degeneracy with the 
peculiar eg transfer amplitudes is important and to be 
a key point . in explaining the various types of AFM 
ordering .E3E2I Until now most of the theoretical studies of 
CO state were done in the framework of the one-otbital 
model ignoring the double degeneracy of eg orbitals.ESO 
The detailed mean-filed analysis of phase diagram of 
one-orbital DE model in the presence of both on-site 
and inter-site Coulomb terms has been given in Ref. pq. 
It has been shown, that in the vicinity of half-doping 
the double-exchange gain of energy is considerably sup- 
pressed by the inter-site Coulomb interaction that favors 
charge-ordered state. Recently, the CO state within the 
two-orbital model has been investigated by the projec- 



tion perturbation techniques combined with the coher- 
ent state formahsm and by Monte Carlo simulations in 
Refs. |l^ and |l^, respectively. In Ref. |l6| the origin of CO 
has been attributed to the effective particle-hole interac- 
tion and CO state with C-type of the spin ordering at 
X — 0.5 has been obtained. The Authors of Ref. |l^ have 
shown that the experimentally observed charge, spin, and 
orbital ordering could be stabilized due to Jahn- Teller 
phonons. 

In the present paper we investigate the role of the or- 
bital degeneracy in the CO state based on the two-orbital 
DE model including the intersite Coulomb interaction. 
We adopt the mean-field (MF) approximation to derive 
the ground state phase diagram in the two-orbital model, 
and compare it to that of the corresponding one-orbital 
model. We argue that the orbital degeneracy together 
with the peculiar e^ transfer amplitude has a drastic ef- 
fect on the phase diagram and is important in obtaining 
the realistic magnetic/chargc/orbital ordering observed 
in half-doped manganites. The paper is organized as 
follows: In the next section the model Hamiltonian is 
presented and the mean-field scheme is formulated. The 
ground state phase diagrams of the one- and two-orbital 
models are derived and compared in Sec. III. Sec. IV sum- 
marizes our main results. In the Appendix the canonical 
transformation diagonalizing the MF Hamiltonian and 
the resulted band structure of various magnetic phases 
are presented. 



II. MODEL AND FORMULATION 

We start with the two orbital ferromagnetic Kondo 
lattice model supplemented by the intersite Coulomb re- 
pulsion 

+ J^SiSj + V^^niftj - ^^nj. (1) 

The first term of Eq.Q describes an electron hopping 
between the two Cg orbitals of the nearest neighbor (NN) 
Mn-ions. The orbitals d'T,^-2_^i and d^i-yi correspond to 
a(/?) — 1 and 2, respectively. Due to the shape of the eg 
orbitals, their hybridization is different in the three cu- 
bic directions that leads to direction dependent hopping 
with the anisotropic transfer matrix elements i" given 



by 



is given by J, V represents the inter-site Coulomb repul- 
sion of Cg electrons, n^ is the particle number operator 
and /i is the chemical potential. The effect of the on-site 
Coulomb interaction that is not included in our model 
Hamiltonian will be discussed later. 

We study the Hamiltonian (|^) within the MF ap- 
proximation, which is set up by introducing the order 
parameter for static charge-density wave of the form 
(jii) = n + (5nexp(iQRi), with n being the electron den- 
sity and Q — (tt, tt, tt). Further, we treat localized spin 
subsystem classically and assume a strong Hund's cou- 
pling 3h ^ ztj S. In this limit one may take the local 
spin quantization axis parallel to i2g-spins and in the ro- 
tated bases retain only " spin- up" components of the mo- 
bile electrons. Then the transfer integral between the NN 
sites is modified through relative angle of the i2g-spins at 

the i and j sites as t°^ = t^^ cos(6'y /2), where % is the 
relative angle of the t2g-spins. We consider the following 
magnetic phases that competes: i) Ferromagnetic config- 
uration (F-type spin ordering) with Q^y = ^z = ( dxy 
and Qz are the angels between the neighboring spin in xy- 
plane and z-direction, respectively, ii) Layer-type anti- 
ferromagnetic configuration (A-type spin ordering) - the 
local spins are parallel in the planes and antiferromag- 
netically aligned between the neighboring planes, that 
corresponds to Qxy = and Qz = ■"'■ hi) Chain-type anti- 
ferromagnetic configuration (C-type spin ordering) - the 
local spins are parallel in the straight chains and antifer- 
romagnetically coupled between the chains - d^y = tt and 
dz — Q. iv) Neel-type antiferromagnetic configuration 
(G-type spin ordering) with all spins being antiparallel 
— dxy ~ Oz = t:. v) CE-type spin ordering with zig-zag 
ferromagnetic chains coupled antiferromagnetically. 
As a result we come to the following MF Hamiltonian: 



{ij),a,l3 
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where A — zV6n, z — 6 for 3-dimensional cubic lattice, 
and d is dimensionality of the magnetic order ( d = 0, 1, 2, 
and 3, respectively for G-, C-, A-, and F- type spin or- 
dering). In Eq.(|^) the zero of the energy is chosen in 
such a manner that the SE energy vanishes in the FM 
state. To obtain phase diagram we need to compare the 
free energies of all possible magnetic configurations. 



III. PHASE DIAGRAM 



A. One-orbital model 



The second term in Eq.(ni) describes the Hund's coupling 
between the spins of localized t2g- electrons Si and the 
itinerant eg electrons with spin ai. The superexchange 
(SE) interaction of localized spins between the NN sites 



In order to incorporate the role of orbital degeneracy, 
first we consider the one orbital model ignoring the dou- 
ble degeneracy of Cg orbitals. Retaining only the one 
orbital per Mn-ion and assuming the isotropic transfer 



amplitude, the electronic part of the MF Hamiltonian in 
k-space is written as: 

^or' - - E(*~k + /^)4rfk - A E 4dk+Q , (4) 



with 



iu = 2t^coski, (ki,k2,k3) = {ka:,ky,k^). (5) 

The above Hamiltonian (H) is easily diagonalized by the 
following canonical transformation: 

dk = "kCi^k + ^kC2,k> ^k+Q = -"^kCi^k + "kCa^k (6) 

with 
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ek = Jil + A2 



(7) 



In terms of the c-operators, the one particle Hamiltonian 
reads 



^ = E(-l)"^(l^)4,aCk,a, «=1,2 



(8) 



At half-filling the chemical potential lies inside the gap 
(/i = 0) and recalling that A = zVSn we receive a self- 
consistent equation for the order parameter 



zV ^ tanh/?ek/2 
2N ^ e^ 



9.N /-^ 



(9) 
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FIG. 1. The zero temperature charge order parameter 5n 
as a function of t/V for the one-orbital model and for a dif- 
ferent (F-, A-, and C-type) spin ordered states. Here t is the 
hopping amplitude and V is the intersite Coulomb repulsion. 



In Fig.|] the overall behavior of the order parameter 
dn as a function of t/V is presented for various mag- 
netic configuration. Since the wave vector summation in 
the right hand side of Eq. (^ diverges in the limit A ^ 
there exist a nontrivial solution even dA, V —^ and hence 
a transition from homogeneous to CO state is continu- 
ous. We also note that Sn diminishes exponentially with 
increasing the bandwidth (see Fig.|l|) indicating that the 
transition between the homogeneous and the CO state is 
a result of the competition between the kinetic and the 
electrostatic energy. 

By comparing the free energies of different magnetic 
configuration we obtain the phase diagram as shown in 
Fig.g. At small V, with increasing J the system, start- 
ing from the F-CO phase, first enters to the C-CO phase 
and then to the G-CO state. Since the gain in the mag- 
netic energy when the system moves from A- to C-phase 
is larger then the gain in the kinetic energy in C to A 
transition the A-CO phase is absent in this part of Phase 
diagram. With increasing oi V at V ~ 0.5t the CO gap 
in C-CO phase overcomes that one in A-CO phase that 
results in opening of small window of A-CO phase in the 
phase digram. We also note that with increasing of V 
the SE coupling needed to stabilize the AFM configu- 
ration decreases since the bandwidth effect is overshad- 
owed when the gap becomes larger. The different mag- 
netic structures in Fig.0 are separated by the first-order 
boundaries. There is a jump in the charge order param- 
eter across the phase boundaries, since the value of order 
parameter dn depends on the effective bandwidth and 
hence on the underlying magnetic structure, as can be 
seen from Fig. 2. 




FIG. 2. Phase digram of the one-orbital model in JS^ /t 
and V/t parameter plane. Here A, C, and G denotes the 
A-, C-, and G-type antiferromagnetic ordering, respectively, 
F denotes the ferromagnetic phase and CO is charge ordering. 
The solid lines stand for the first-order phase boundaries. 



B. Two orbital model 



To describe the effect of orbital degeneracy, we con- 
sider the MF Hamiltonian (H) with anisotropic hopping 



amplitude. In the momentum space the electronic part 
of the Hamiltonian reads as: 
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with 
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(11) 



i COS Uxy-, T^z 



t COS 6z ■ 



ordering, there is no difference between the one and two 
orbital models in this sense. The existence of the addi- 
tional orbital only introduce the localized level which is 
empty at the filling we consider. 

The above Hamiltonian (10) describes the F-, A-, C-, 
and G-type spin ordering on an equal footing. How- 
ever the presence of additional orbital degree of freedom, 
with the peculiar anisotropic transfer amplitudes t"^ _ 
[see Eq.(||)] results in the anisotropic DE interactionHl 
and may lead to the stabilization of CE spin ordering. 
Let us consider one zig-zag with two ferromagnetic bonds 
alternated in x and y directions (Fig.0). The corner and 
middle sites are denoted by a{a) and 5(6), respectively, 
and the unit cell is given by four nonequivalent atoms. 
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FIG. 3. Order parameter Sn versus t/V for the two-orbital 
model and for a different (F-, A-, and CE-type) spin ordered 
states. 

The diagonalization of this Hamiltonian ( see Ap- 
pendix) leads to the four band model. In the case of F- 
and A-type spin ordering and at the filling correspond- 
ing to one electron per two Mn-ions (half-doped case) the 
gap is not opened at the Fermi surface, and the chemi- 
cal potential moves down with the lower two bands. We 
solve the gap equation self-consistently with one for the 
chemical potential. As it seen in Fig.^, the transition 
to the charge ordered state is not continuous and there 
exists a critical value Vc above which the ordered state 
is favorable (V/ ~ 0.72i and V/^ ~ 0.58i for F- and A- 
type spin ordering, respectively). As for the C-type spin 
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FIG. 4. The one isolated zig-zag shaped ferromagnetic 
chain. The real CE-type spin structure is given by this 
quasi-one-dimensional ferromagnets coupled antiferromagnet- 
icaly. White and black circles denote Mn**^ and Mn^^ ions 
respectively. The dashed line shows the unit cell and a' is the 
lattice parameter along the zig-zag 

For further discussions it is convenient to adopt the foL 
lowing bases of the eg orbitals at nonequivalent sites:E£l 
|1> = d3^2_^2, |2) =_dy2„^2, and |1) = d3y2_^2 12) = 
dx2_z2 on a, b and a, b sites, respectively [see, Fig.^. In 
the new bases the transfer matrix elements are given by 



t' 



Q/3 



t 



1 




Q/3 



t 



-1/2 
\/3/2 



(12) 



between a-b ( a-b) and b-a ( b-a) NN sites, respectively. 
As a result the zig-zag chain is modeled as a dimerized 
one with the alternating hopping amplitude and can be 
described by the following Hamiltonian 



TirCE _ 



E { [^I'^i^Lbtp + alabtfi} + t'i'^iblaatp + b\^a^+ii}} + Ha 



i,a.f3 



ifi + A){al^a^a + 4aO■^<^} - (a* - A) {bj^b.^ + ^L^*"}} 



(13) 



where i runs along the zig-zag and denotes the number 
of unit cell. Diagonalization of the above Hamiltonian 
( see Appendix) leads to the complicated band struc- 
ture consisting of bonding and antibonding bands, Ea.b — 
±-\/A2 + ^2(2 — cos(fc/2)), and nonbonding states E± = 



±A. Due to the topology of the zig-zag structure, only 
the directional d33,2_^2(d3y2_^2) orbitals at a{a) sites give 
the input in the low energy bonding state [see Appendix 
for a details]. While for the b(b) sites both two orthog- 
onal orbitals do contribute. Therefore, the orbital de- 



generacy is removed on the middle site sublattice and 
the carriers on this sublattice will occupy the direc- 
tional orbitals, leading to the polarized orbital state with 
d^^2 _^2 / d^y2 _j.2 orbital ordering. We also emphasize, 
that at half-doping the bonding band is full and the sys- 
tem is a band insulator even in the absence of the charge 
ordering. The onset of the charge ordering renormalizes 
the gap to higher value. The behavior of the charge order 
parameter is depicted in Fig. 3. The transition to charge 
ordered state takes place at V^^ ~ 0.44i, that is lower 
then that one for A- and F- type spin ordering. A smaller 
value of intersite Coulomb repulsion is needed to intro- 
duce the charge ordered state in the state with the lower 
dimension of ferromagnetic component. 
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FIG. 5. The phase diagram of the two-orbital model {JS jt 
vs V/i4) . The dashed lines stand for the second-order bound- 
aries and separate the uniform state from the charge ordered 
one in the corresponding magnetic phases. 



The ground state phase diagram of the two-orbital 
model is given in Fig. 5. The phases with different mag- 
netic structures are separated by the first-order bound- 
aries (soled lines in Fig. 5). The transition from the 
charge-disordered to the charge-ordered version of the 
given magnetic structure is continuous, there is no jump 
in the charge order parameter during the transition which 
does not change the symmetry of underlying magnetic 
structure ( the corresponding second-order boundaries 
are denoted by the dashed lines in Fig. 5). At small 
value of V with increasing of SE coupling the system 
goes through four magnetic phases F, A, CE, and G con- 
sequently. Due to the additional orbital degree of free- 
dom, unlike to the one orbital model, there exist a finite 
phase space of A-type spin ordering in the small V re- 
gion V < V^. However, for V > V/^, when the charge 
ordering is introduced in the A phase, and hence the 
bandwidth effect is suppressed, CE-CO phase wins and 
the A-CO phase is never realized in the phase diagram 
of the present model. 



IV. CONCLUSION 



In summary, we have investigated the ground state 
phase diagram of the one- and two-orbital extended dou- 
ble exchange models within the mean filed approximation 
at half-doping. In the case of the one-orbital model the 
MF theory predicts the continuous phase transition to 
the charge ordered state due to the intersite Coulomb 
interaction V. In the two-orbital model the character of 
the transition is changed by the introducing the nonzero 
critical value of V. While the transition to charge ordered 
state with given magnetic structure is continuous, there 
is jump in of charge order parameter across the phase 
boundaries between the states with different symmetry of 
magnetic structure. Depending on the intersite Coulomb 
interaction V and and superexchange coupling J differ- 
ent types of spin ordering (F-, A-, G-type) accompanied 
by the charge ordering may take place in the ground state 
of the one-orbital model. The presence of orbital degen- 
eracy with the peculiar anisotropic eg transfer amplitude 
introduces the new magnetic state, CE-type spin order- 
ing, in the phase diagram of the two-orbital model. The 
C-type spin ordering is never achieved within this model 
due to its instability against the effective " dimerization" 
and formation of the zig-zag ferromagnetic order. The 
alternation of the ferromagnetic bonds in x/y directions 
leads to the alternation of the hopping amplitude. As 
a result the bare band is splited into bonding and anti- 
bonding states and the "dimerization" gap opens on the 
Fermi surface at half-doping. The CE-type spin ordered 
state is accompanied with d^^'^-r^ /d'iy2_j.2 orbital order- 
ing, originated from the topology of the zig-zag structure. 
The CE-type spin/charge state also wins the CO state 
with A-type AFM ordering and there exists the phase 
boundary between the CE spin/charge ordered state and 
charge disordered (CD) A-type state [see Fig. 5]. The 
experimentally detected competition between A-CD and 
CE-CO states in the half-doped manganites could indi- 
cate that the parameters of the system are close to this 
phase boundary. Therefore the small change of the band- 
width may drive the system from one to other state. That 
is also suggested by the huge, oxygen isotope effect ob- 
served in La-Pr compounds .E3'l2I A substitution of O^^ 
by the isotope O^^ narrows the carrier bandwidth due to 
the polaronic effect and as a result the electrostatic en- 
ergy might overcome the kinetic energy and the charge 
ordered insulating state might be established. 

We would like to point out that there exist other phys- 
ical factors and ingredients not included in present treat- 
ment that may stabilize the particular state and modify 
the phase diagram (quantum nature of the spins, cou- 
pling to the Jan- Teller phonons, and on site Coulomb 
interaction) . As it was recently discussed in Ref. |l] the 
interorbital on-site Coulomb interaction U' could be re- 
sponsible for the experimentally observed charge ordered 
state. In the CE structure the orbital degeneracy is re- 
moved on the sublattice composed by the middle sites 



and the low energy state corresponds to the directional 
orbitals (see Appendix of the present paper). If charge 
carriers occupy the corner sites there will be a positive 
contribution from the U' term to the system energy. 
While, if only middle site sublattice is occupied this pos- 
itive contribution disappears since the onsite Coulomb 
term has the zero matrix element in this case. Therefore, 
U' term in CE-type spin ordered state acts as a source of 
the charge ordering, but in a rather deferent way then the 
intersite Coulomb repulsion. Since the middle sites are 
stacked in ^-direction U' induces experimentally observed 
(tt, tt, 0) charge ordering instead of the Wigner crystal 
type CO favored by the intersite Coulomb interaction 
considered here. Nevertheless, we believe that the main 
effect of the orbital degeneracy on the interplay of differ- 
ent type of magnetic ordering in charge-ordered state is 
at least qualitatively captured by the present treatment 
indicating the crucial importance of the electronic state 
degeneracy on the phase diagram. 
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One finds the effective Hamiltonian of the form 



kA k 



(A3) 



Further we perform the transformation similar to (|6|) and 
introduce the new sets of fermionic operators as 

ftkl = Wkl^kl + '^klrfk+Q27 ak2 = '"klrfk+Q2 ^ ^kl^kl; 
%3 = "k2rfk2 + 'i'k2rfk+Ql) ^^kA = '"k2rfk+Ql ^ ^k2rfk2: 

where 
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APPENDIX: DIAGONALIZATION IN 
MOMENTUM SPACE 



1. Two orbital model 



Diagonalization of the Hamiltonian (nffl) can be done 
by two subsequent canonical transformation. First, we 
diagonalize the free part of the Hamiltonian by introduc- 
ing the new fermionic operators 

{ dki = Ukdki - Vkdk2, dk2 = -Vkdki + Ukdk2 (Al) 



with 



"ki 



V2 



Ski 



Vki 



Ek^ = ^|eh + A\ 



f .N» sgn(4^) 
i = 1,2. 



£ki 



As a result we get the following four band Hamiltonian 



rr2oib 



2J |-Eki[aj^iaki - ak2ak2] 

k 

4 

£^k2[ak3"k3 - ak4"k4]| " M X! "ki«ki- (A4) 



k,i=l 



2. CE-structure 



First, we rewrite the linearized MF Hamiltonian ( [Ti^ ) 
in the momentum space: 



TirCE _ 



ir{&L«k/5 + eni,„ak/5} + i?. 



kj;t 

ka"^ 



Z! { \tfWkoPkii + aL^k/3} 

k,a,/3 

{^i -I- A){a[^aka + aLoka} - {p - A) {^L^ka - ^L^ka}| 



(5) 



where —Tr<k<TT and the lattice constant is set to 
be unity (a' = 1 see Fig.^). The above Hamiltonian can 
be simplified by transforming to new fermion operators 

\ ska 1 Ska i ^ka ; ^ka / 5 



The operators {rjk, Vk}, are obtained from 6ka, ^ka by the 
same transformation as in Eq.(|6|). One finds the effective 
Hamiltonian of the form H^^ = H + H, where 



flko 



ska t ska 

V2 ■ 



a-kc 



= e*/2 



ska I ska 

V2 ■ 



(6) 



H= -T. {tf^Lvw3 + tfe^/'vU^, + +H.c] 
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H ^ H 






^)^L?ka + (M-A)?7L'?k«}, (7) 
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ap 



(8) 



Further, wc consider only the first part given by (M), 
generalization of the diagonalization procedure on H is 
straightforward. With the help of the explicit expression 
of the hopping amplitude matrix we first transform oper- 
ators 77ki.2 to ?7ki,2 by the fermionic u — v transformation 
with 



,ik/2 



Mk 



/2 



^v 



I'k 



\/3e*/V2 



fik 



(9) 



where fik = ■\/2 — cosk/2. In terms of the new operators 
the effective Hamiltonian is written as 



H = - 



fk{eli77ki + H.c.} 
^ {(m + A)fi„Ck. + (m - A)r)t^77ka} . 



k 



(10) 



ka 



where £k = if^k- The transformed Hamiltonian is al- 
ready diagonal in the subspace given by the effective or- 
bital index a = 2. However it mixes the fermionic fields 
Cka and 77kQ at a = 1. The hybridization part of the 
Hamiltonian (|lO|) can be diagonalized following the same 
rout as in Eqs. (Al A2). As a result we come to the 
following diagonal form: 

J? = EK[/3li/3ki-/3lA2] 

k 

4 



k,i=l 



H = H[I3^ I3,8^£] 



(11) 



where 



fk = VA2-^t2(2-cosfc/2), 
4 = VA2-Ft2(2-f cosfc/2), 



(12) 

and /3(^)k3 = '7(f])k2, /3(;5)k4 = C(C)k2, and /9(^)ki,2 are 
given by the linear combination of C('?)ki and ?7(?7)ki- 

Let us now briefly comment the obtained band struc- 
ture. It is given by the bonding and antibonding bands 
=F£(£)k ( p2[ ) and nonbonding bands ±A in between. The 
nondirectional orbital from the middle sites {dy2_z2 and 
d^2_^2 orbitals on a and a type sites, respectively) is 
completely decoupled from the other states giving rise 
to the nonbolding band with energy —A corresponding 
to /3(/3)k4 in Eq.(|l^). The other nonbonding band with 
energy A corresponding to the states /?(/3)k3 is given 
by the linear combination of the degenerate orbitals on 
the corner sites. The state orthogonal to this nonboning 
state hybridizes with the directional orbital on the mid- 
dle sites (d3a;2_j.2 and djy2_j.2 orbitals on a and d type 
sites, respectively) leading to the bonding and antibond- 
ing bands. As follows from the above discussion, the par- 
ticular geometry of the zig-zag structure not only leads to 



the opening of "dimerization"-like gap in the spectrum, 
but also removes the orbital degeneracy on the middle 
sites. The energy of the directional orbital is lowered 
due to the hybridization and hence delocalization, while 
its orthogonal orbital remains local. 
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